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DECUS PROGRAM LIBRARY WRITE-UP DECUS No. FOCAL8-69 

TITLE: ANALYSIS OF VARIANCE 

ABSTRACT: This program is designed to solve the analysis of 
variance problem for the two-factor completely randomized 
design, and to table. the results of the analysis in a form 
acceptable for publication in many scientific journals. Both 
the input and output formats are designed for simplicity and 
ease of operation. 

An alternate form of the program makes possible the evaluation 
of either one-factor or two-factor designs. 

OPERATION: Both forms of the program must be loaded without 
the extended functions. The following equations are used: 

C (a calculated correction term) : ¥\7rET~™ 



A (row variable) 
B (column variable) 






AB (interaction) : ^fzl±^Li _ £ - A~*& 

T (total variance) : j? ££ y 2 — n 

S/AB (within cells variance) : T~C~A"*6~AB 

The program is initiated by typing "GO" followed by a 
carriage return. The carriage return is the terminal character 
for all questions in the program. The pro gram will then ask 
"N" which should be answered by the number of entries in each 
cell (the program is designed only for equal values in each 
cell, hence "N" is asked only one time). This is followed by 
questions of "No. of Rows" and' No. of Columns"which should be 
answered appropriately. Input is then made into the samples 
which are identified by a signed pair number of the following 
form: (R,C) where R=the row number of the sample, and 
C=the column number of the sample . 

By responding with a 1 to either "No. of Rows" or "No. of 
Columns" in Format Two, the program performs a one factor analysis. 

SPECIAL OPERATION: The following modifications to FOCAL 
(DEC-08-AJAD-PB) should be made in order for the input and 
output formats to retain their neat appearance. 



1261/0200 
1262/0200 



: Eliminates printing of lf = fl 
: Eliminates printing of ":" 



RESTRICTIONS: The maximum size of input into Format One is 
approximately 27 cells of data. The use of Format Two reduces 
this to approximately 20 cells. The number of entries within 
each cell (N) is for all practical purposes, unlimited. 

The output format is set at %7.02, i.e., seven significant 
digits with two places to the right of the decimal point. There 
may be occasions when this will cause some loss of accuracy, 
and the format may be changed to floating point by modification 
of the appropriate statements in command section 3.0 in either 
program. 



ANALYSIS OF VARIANCE- FORMAT ONE 

C-FOCAL , 8/68 

0\.\0 ERASE 

0\.20 A "N "P;A "NO. OF ROWS "H 

jell. 21 A "NO. OF COLUMNS 'V,!'. 

0].30 F N=1,H;D0 4 

02A0 F B=1,H;D0 7 

02.20 F BB=1,V;SSS=BB-V;D0 9 

02.25 F B=1,H;S J1=E(B)+J1 

J02.3J0 S C=Jlt2/(H*V*P) 

02.35 F B=1,H;S J>E(B)f2+JJ 

,02.40 S Vl= [JJ/(V*P)]-C 

J02.45 F BB=1,V;S J2=EE(BB)T2+J2 

02.50 S V2=LJ2/(H*P)]-C 

J02.65 F U=1,H*V;S L=A(U)f2+L 

J02.7J0 S V3=[L/P>C-V1-V2 

02.75 F >1,R;S J3=G(J)+J3 

JEI2.81 S T=J3-C 

J02.9J0 S ER=T-V1-V2-V3 

03. \0 T .'.'.'. , / "ANALYSIS OF VARIANCE TABLE", 1 1 i 

03.20 T "SOURCE SS DF MS F",il 

J03.3J0 T %7.J02 TOTAL '7, (H*V*P)-l,li 

.03.40 T "ROWS ",V1, H-l, V1/(H-1), (Vl/(H-l))/<ER/CH*V*(P-lfl> / , .i 

03.50 T "COLUMNS ",V2, V-l, V2/(V-1), (V2/(V-l))/<ER/fH*V*(P-l) = ]> / , .l 

J03.6J0 T "R X C ",V3, (H-1)*(V-1), V3/<(H-1)*(V-1)> 

03.65 T (V3/<(H-1 )*(V-1)>)/<ER^H*VMP-1 )]>, '. 1 

03.70 T "ERROR VAR ",ER,£h*V*(P-1)J, ER/<H*V*(P-l)>,Ii 

03.80 QUIT 

J04. 1J0 F 0=1 , V;S R=R+1 ; DO 5 

J05 . 1J0 T %1 , "SAMPLE ( " N, ", ",0, "), \ 

05.20 F Q=1,P;D0 6 

06.\0 A X(R);S A(R)=X(R)+A(R);S G(R)=X(R)f2+G(R) 

07. \0 F D=1,V;S S=S+1;S E(B)=A(S)+E(B) 

09. ]0 F DD=1 / H;SSS=SS+V;SEE(BB)=A(SS)+EE(BB) 



* 



N 8 

NO. OF ROWS 2 

NO. OF COLUMNS 3 

SAMPLE ( I, 1) 

7 

33 

26 

27 

21 

6 

14 

19 

SAMPLE (1, 2) 

6 

11 

11 

18 

14 

18 

19 

14 

SAMPLE (1, 3) 

9 

12 

6 

24 

7 

1 

1* 

SAMPLE ( 2, 1) 

42 

25 

8 

28 

3j* 

22 

17 

32 

SAMPLE ( 2, 2) 

28 

6 

1 

15 

9 

15 

2 

37 



SAMPLE ( 2, 3) 

13 

18 

23 

1 

3 

4 

6 

2 



ANALYSIS OF VARIANCE TABLE 



SOURCE 


SS 


DF 


MS 


F 


TOTAL 


4987.91 


47.0 






ROWS 


40.33 


i.jejtf 


40.33 


0.49 


COLUMNS 


1387.J04 


2.0 


693.52 


8.49 


RXC 


127.54 


2.0 


63.77 


0.7B 


ERROR VAR 


3433 .0 


42.0 


81.74 





ANALYSIS OF VARIANCE- FORMAT TWO 

C- FOCAL , 8/68 

01.10 ERASE 

161.20 A "N "P;A "NO. OF ROWS "H 

J01.21 A "NO. OF COLUMNS "V,ll 

01.30 F N=1,H;D0 4 

j02.1jef F B=1,H;D0 7 

02.20 F BB=1,V;S SS=BB-V;DO 9 

J02.25 F B=1,H;S J1=E(B)+J1 

02.30 S C=J1T2/(H*V*P) 

J02.35 F B=1,H;SJ>E(B)T2+JJ 

J02.40 S V1=[JJ/(V*P)>C 

J02.45 F B&=1,V;S J2=EE(BB)T2+J2 

02.50 S V2=[J2/(H*P)>C 

J02.65 F U=1,H*V;S L=A(U)T2+L 

02.70 S V3=[L/PJ-C-V1-V2 

02.75 F J=1,R;S J3=G(J)+J3 

J02.81 S T=J3-C 

02.^0 S ER=T-V1-V2-V3 

0Z.10 T 1 i 1 I, "ANALYSIS OF VARIANCE TABLE", 1 1 1 

03.20 T "SOURCE SS DF MS F",il 

Jd3.3,0 T %7.02 "TOTAL "T, (H*V*P)-1,11 

,03.35 IF (H-l) 3.4,3.5,3.4 

,03.40 T "ROWS ",V1, H-l, V1/(H-1), (V1/(H-1))/<ER/£h*V*(P-1 )>,'.'. 

03.45 IF (V-l) 3.5,3.7,3.5 

J03.5J0 T "COLUMNS "V2, V-1, V2/(V-1), (vy(V-l))/<ER/CH*V*(P-l)]>,ll 

J03.55 IF ((H-1)*(V-1)) 3.6,3.7,3.6 

J03.6J0 T "RXC ",V3, (H-1)*(V-1), V3/<(H-1)*(V-1)> 

J03.65 T (V3/<(H-l)*(V-l)>)/<ER/[;H*V*(P-lO>, , .l 

03.70 T "ERRORVAR ",ER, QH*V*(P-lJ], ER/<H*V*(P-1)>, 11 

J03.8J0 QUIT 

J04.1J0 F 0=1,V;SR=R+1;D0 5 

05.10 T %1 , "SAMPLE ("N, ", ", O, ") ", 1 

05.20 F Q=1,P;D0 6 

06.10 A X(R);S A(R)=X(R)+A(R);S G(R)=X(R)T2+G(R) 

07.10 F D=1,V;SS=S+1;SE(B)=A(S)+E(B) 

J09.1J0 F DD=1,H;S SS=SS+V;S EE(BB)=A(SS)+EE(BB) 



*GO 

N 8 

NO. OF ROWS 1 

NO. OF COLUMNS 5 



SAMPLE (1,1) 

5.8 

5.1 

5.7 

5.9 

5.6 

5.4 

5.3 

5.2 

SAMPLE (1, 2) 

6.0 

6.1 

6.6 

6.5 

5.9 

5.9 

6.4 

6.3 

SAMPLE (1,3) 

6.3 

5.5 

5.7 

6.0 

6.1 

6.2 

5.8 

5.6 

SAMPLE (1, 4) 

6.4 

6.4 

6.5 

6.1 

6.6 

5.9 

6.7 

6.0 

SAMPLE (1, 5) 

5.7 

5.9 

6.5 

6.3 

6.2 

6.4 

6.0 

6.3 



ANALYSIS OF VARIANCE TABLE 

SOURCE SS 

TOTAL 6.33 

COLUMNS 3.48 4.00 0.87 10.72 

ERROR VAR 2.84 



DF 


MS 


29.00 




4.00 


0.87 


35.00 


0.08 



